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Abstract: We compute next-to-leading order (NLO) corrections in the e-regime of Wil- 
son (WChPT) and Staggered Chiral Perturbation Theory (SChPT). A difference between 
the two is that in WChPT already at NLO, that is at 0(e 2 ), new low energy constants 
O ; (LECs) contribute, whereas in SChPT they only enter at C(e 4 ). We first determine the 

NLO corrections in WChPT for SU(2), and for U(Nf) at fixed index. This implies finite- 
volume corrections to the phase boundary between the Aoki phase and the Sharpe-Singleton 
scenario via corrections to the mean field potential. We also compute NLO corrections to 
the two-point function in the scalar and pseudo-scalar sector in WChPT. Turning to SChPT 
we determine the NLO corrections to the LECs and their effect on the taste splitting. Here 
the NLO partition function can be written as the leading order one with renormalized cou- 
plings, thus preserving the equivalence to staggered chiral random matrix theory at NLO 
for any number of flavors Nf. In WChPT this relation only appears to hold for SU(2). 
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1 Introduction 

In the study of the low energy dynamics of QCD it has become of considerable interest 
to extend chiral perturbation theory by including effects of the lattice spacing a. Indeed 
in this way one can control how the ultraviolet cut-off influences every numerical simu- 
lations of lattice QCD. For what concerns the Wilson fermion formulation, in a series of 
works [1, 2, 3, 4, 5] it was shown how the continuum chiral Lagrangian gets modified from 
discretization effects, so-called Wilson Chiral Perturbation Theory (WChPT). For the stag- 
gered version of QCD the analogous extension has been done by Lee and Sharpe [6] for 
one staggered flavor and then generalized by Aubin and Bernard [7], so-called Staggered 
Chiral Perturbation Theory (SChPT). 

Quite recently there has been a lot of activity in the study of the e-regime of WChPT and 
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SChPT. The discretization effects compete with the quark mass for the explicit breaking of 
chiral symmetry. Despite the fact that this regime originally introduced in the continuum 
in [8] is unphysical, which is due to the fact that Compton wave-length is larger than the 
size of the box, it can be extremely useful in the determination of the low energy effective 
constants (LECs) and in the study of the Dirac operator spectrum and its dependence on 
the topology of the gauge fields [9, 10, 11, 12, 13, 14, 15, 16]. Due to the domination of the 
zero-mode group integral many computations can be performed analytically. First lattice 
simulations have already confirmed the quenched predictions [17, 18, 19]. 
Another intriguing characteristic of this regime is that at leading order (LO) in the e- 
expansion it is equivalent to a chiral Random Matrix Theory (ChRMT). More precisely, in 
the continuum theory a rigorous proof of this equivalence has been given in [20, 21] for all 
Dirac operator eigenvalue correlation functions. The extension to theories at finite lattice 
spacing have been derived subsequently: in [10, 11, 12, 13, 14, 15, 16] the authors showed 
that an analogous equivalence holds between the LO of WChPT and a Wilson Chiral Ran- 
dom Matrix Theory (WChRMT). In [22] this equivalence between the zero-mode sector of 
SChPT and the corresponding Staggered Chiral Random Matrix Theory (SChRMT) was 
established. 

From the point of view of ChPT when introducing the effects of finite lattice spacing one 
has to consider the symmetry breaking of the continuum theory down to a subgroup. Con- 
sequently more operators will be allowed in the Symanzik effective action, and more terms 
appear in the chiral Lagrangian. In addition to the chiral condensate E and to the pion 
decay constant F some new LECs need to be introduced to characterize the strength of 
these terms. In WChPT in general 3 new LECs are introduced at LO and labeled by 
Wq, Wj, Wg, whereas for SChPT 6 new LECs need to be introduced, C\, C2, C3, C4, C5, Cq. 
The absolute sign of individual LECs and of their combinations has been subject of intense 
recent discussions in WChPT [11, 23, 24, 25]. It remains to be seen if such arguments carry 
over to NLO. 

In the e-regime there are different ways of introducing the relative strength between the 
mass terms and the lattice spacing. In this paper we will use the so-called Aoki or large 
cut-off effect (LCE) counting where m ~ a 2 AQ CD or a ~ 0(e 2 ) [26]. Thus the discretiza- 
tion terms are of the same order as the mass terms and compete for the breaking of chiral 
symmetry already at LO. Considering the Aoki regime means that the LO integral over the 
zero-modes is modified with respect to the continuum theory and makes the analytic com- 
putations more involved. However, one could also consider the Generic Small Mass (GSM) 
counting [27, 28] where the discretization errors enter only at next-to-next-to-leading order 
(NNLO) with respect to the continuum Lagrangian, since the counting is m ~ aAq CD 
or a ~ C(e 4 ). Finally an intermediate regime is known as the GSM* counting, where 
a ~ C(e 3 ) and the discretization errors are at next-to-leading order (NLO) with respect to 
the continuum [29, 30]. 

The main question addressed in this paper is to extend WChPT and SChPT to order C(e 2 ), 
for both the partition functions, and for the scalar and pseudoscalar current correlators in 
the Wilson case. For the spectral density of the Dirac operator a NLO order calculation 
has already been done, however not in the e- but in the p-regime [31]. In contrast to the 
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e-regime in the continuum, sectors of fixed topological charge v are no longer well defined 
at finite-lattice spacing. They have to be replaced by the index of the Dirac operator, and 
we refer to [11] for a detailed discussion and references. In the continuum it was found 
that the NLO partition function could be written as the LO order one with renormalized 
couplings that include the order 0(e 2 ) effects, both for £ and for F when including a iso- 
spin chemical potential [32, 33, 34]. This implied that the partition function from ChPT 
and from ChRMT agree up to NLO. Only at NNLO non- universal effects were found in 
[35]. We will find that for WChPT only for SU(2) the NLO contributions are absorbed 
into the two effective couplings relevant in that case, whereas for U(Nf) at fixed index we 
have to take extra derivatives of the LO partition function. In contrast for SChPT the 
NLO contributions can be absorbed into renormalized couplings for any Nj, due to the 
U{\) remnant of the continuum chiral symmetry. 

Other information can be extracted from the NLO partition function. For SU(2) one can 
see that the finite-volume corrections can drive the system in or out of the Aoki phase, 
compared to LO. Depending on the coefficients of the mean field potential the theory can 
stay in the Aoki phase where two pions are massless as a consequence of the breaking of 
the flavor symmetry, or in the Sharpe-Singleton scenario where a first order transition is 
present. The finite-volume effects can modify the boundary of these two regions. 
For what concerns staggered fermions the effective LECs we compute can lead to the fol- 
lowing prediction. Since from the tree level Lagrangian one can see how the taste symmetry 
is broken, the new renormalized LECs allow to quantify how the finite-volume corrections 
can modify the taste symmetry violation. 

The outline of the paper is as follows. In section 2 we study WChPT at NLO starting with 
SU(2) and including the effect on the phase boundary in subsection 2.2, and then turn to 
U(Nf) in 2.3. The two-point functions are given in 2.4. In section 3 we repeat our analysis 
for the staggered version, which includes the effect on the taste splittings. Finally in the 
last section 4 our discussion and some considerations regarding possible extensions of this 
work are presented. Several technical details are deferred to the appendices A to E. 

2 Wilson Chiral Perturbation Theory at NLO 
2.1 Introduction 

In this section we consider the e-regime of Wilson Chiral Perturbation theory (WChPT) 
with Nf = 2 degenerate quarks of mass m. As already pointed out in the previous section 
and shown in [11] at LO it is equivalent to a Random Matrix Theory which includes 
order 0{a 2 ) discretization effects (WChRMT). Our aim is to analyze WChPT at the next- 
to-leading order (NLO) in the e-expansion and show that within the Aoki regime the 
partition function at that order can be rewritten as the LO one with renormalized low 
energy constants (LECs). In the continuum a similar relation between the LO and NLO 
partition function holds for every number of flavors as shown e.g. in [32, 33]. 
Let us start by introducing the two-flavor Wilson chiral Lagrangian that at LO in the Aoki 
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regime can be written as 
F 2 



-Clo 



-Tr 



-Tr 



M t f7+ U^M 



+ a 2 C2 



Tr 



(2.1) 



In addition to the continuum Gasser-Leutwyler terms [36, 37] there is an additional order 
0(a 2 ) contribution and thus a new low energy effective constant C2 (note that our c 2 = 
Wq + Wg/2 is a short hand notation for the standard terminology which is C2-F 2 /16, and 
likewise for the other NLO LECs). Here as usual F is the pion decay constant, S is the 
chiral condensate and M is the mass matrix that, for two degenerate quarks with mass m, 
reduces to M = ml.2y.2- 111 an the following we will only consider degenerate masses. We 
use the standard parameterization for the Goldstone boson 



U (x) = Uq exp 



(2.2) 



where Uq is the two by two unitary matrix describing the zero-modes nonperturbatively, and 
the Hermitian fields £(x) = = <7ft&>, that belong to the Lie algebra su{2), parameterize 
the propagating modes. In order to derive the Lagrangian (2.1) in the so-called Aoki regime 
one has to use the power counting [26] 



V ~ e 



m ~ 



e 4 , d ~ e, £(x) ~ e, a ~ e 2 



(2.3) 



Other counting schemes can be considered if we want to study the GSM* or GSM regime 
[27, 28] that are in fact defined considering the cut-off effects respectively as NLO and 
NNLO contributions with respect to the continuum terms. 

If we want to go further in the Aoki regime and compute the NLO partition function we 
have to consider in addition to (2.1) the NLO chiral Lagrangian. 



Leading Order O(e ) 


2 2 

m, p , a 


Next-to-Leading Order 0(e 2 ) 


am, ap 2 , o 3 


Next-to-Next-to-Leading Order 0(e 4 ) 


2 2 4 2 2 2 4 

m , mp , p , a m, a p , a 



Table 1. Contributions to the Wilson Chiral Lagrangian in the Aoki regime. 

As we can see from the table 1. that schematically gives us the counting of all the terms 
that contribute up to C(e 4 ) to the chiral Lagrangian, at order 0(e 2 ) the possible terms 
that enter are 0{ap 2 ), O(am) and 0(a 3 ). Following [2, 24, 38, 39] they can be written as 

2 



£nlo = a c Tr 

+a 3 d 1 Tr 



Tr 



u + u ] 



+ a 3 d 2 ( Tr 



+ am C3 ^Tr 
+ 1\3 



(2.4) 



where for SU(2) 4 new and undetermined LECs, namely Co, C3, d\ and a\, need to be intro- 
duced. Here we have explicitly used some special properties of SU(2), see e.g. appendix 
A, compared to the general Nf case. Now using the power counting (2.3) we expand the 
action 

S = d A x (C LO + £nlo) (2-5) 
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up to 0(e 2 ), where we obtain 

S-(o) = I [ d 4 x Tr [d^(x)d4(x)} - \ m Tr [t7 + E/J 



+ a 2 Vc 2 Tr 



<j(0) , q(0) 



(2.6) 

for the O(e ) contribution. Here we have defined the LO part of propagating and zero- 
modes separately. For the 0(e 2 ) terms we get 



5(2) 



12 F 2 



d xTr 



[d^{x),t(x)][d^{x)^{x)] 



+ 



2F 2 



d 4 xTr 



U + c/ f ) £(z 



-2a 2 1| i 



+ 



2a cq 



Tr 



F 2 

+a 3 diFTr 



o + ^ 
Co + Ul 



d 4 xTr 



Co + Ul £(x) 



U + U< 



d 4 x Tr [d^(x) d^(x)] + amc 3 V ( Tr 

3 



Uq + U, 



+ a 6 d 2 V Tr 



U + U, 



(2.7) 



while the 0(e) term vanishes due to f d x£(x) = 0. 

2.2 Partition function for Nf = 2 and Aoki Phase 

The next step is the calculation of the partition function up to 0(e 2 ). The general form of 
the partition function can be rearranged by separating the integration over the zero-modes 
from the integration over the Gaussian fluctuations as 



Z 



SU(2) 



[d H U(x)} e 



,,(0) 



SU{2) 



d H U e-°^o Z^U ) , 



2.i 



with 



Z € (U ) = J [dt{x)] (l - ^ J d'xTr [^(x) 2 ]) e S %~ S , (2.9) 

containing the Jacobian J(£(cc)) up to order C(e 3 ) [8] from the parameterization (2.2), and 
the chiral action S to an unspecified order. Here the invariant Haar measure 

N f 



[d H u( X )] = d H u [de(x) 



1 



d 4 xTr [f (a 



(2.10) 



3F 2 V 

has been divided as the invariant measure over the zero-modes Uq times the fiat measure 
over the fluctuation £. At this point one can expand the function Z^(Uo) up to 0(e 2 ) and 
then perform all the Gaussian integrals using the expression 

1 



[d£(x)] expf-S 1 ^] €(x)ij£(y)ki = 

in terms of the propagator. We easily find that 

3mVZ 



8u$jk - jy- 5ij8 k i ) A(x - y) 



(2.11) 



Z^Uq) = MU + 



4F 2 



A(0) - a 3 diV Tr U + U% 



+ 



4a 2 c 2 V 
F 2 



A(0) - amc 3 V Tr 



Uo + U, 



a 6 d 2 V Tr 



U + U, 



(2.12) 
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Here M is an overall normalization factor that contains all constants that are Uq indepen- 
dent and that drop out in expectations values. In particular this includes the contribution 
from the Jacobian. The propagator A(0) is finite in dimensional regularization and is given 
by A(0) = — fii/V 1 / 2 , with (3\ a numerical coefficient that encodes the geometrical data 
of the box considered. At this point we note that all terms in (2.12) can be reabsorbed 
easily in the LO chiral Lagrangian by re-exponentiating the corrections, with the only the 
exception of the last term. In order to solve this problem one can write this contribution 
as a sum of single and double trace terms using the relation (A. 8) obtained in the appendix 
A. 2 through some group integral identities. Finally the partition function can be written 
as 



-^nlo = M ' d H U exp 

JSU{2) 
_ N 7 / v eff fiffN, 



ms cff y 



-Tr 



U + Ul 



a 2 cf V[Tv 



U + Ul 



with the effective renormalized LECs given by 



SI 



2F 2 



A(0) 



771V V 



2a 2 d x + 32a 2 d 2 - 3— 

C2 



and 



r eff 



C 2 1 



F 2 



A(0) + 



m 



C3 d 2 
S 4c 2 



V 



Here we have defined 



m 



mT,V and a = a V 



(2.13) 



(2.14) 



(2.15) 



(2.16) 



which are of order 0(1). Differently from the continuum limit chiral perturbation theory, 
in which the NLO renormalized LECs can be rewritten only as functions of the number 
of flavors and the geometry of the system, here their expressions involve also some NLO 
LECs. In principle this allows us to extract them from lattice computations through a 
finite-size scaling analysis. Performing the simulations at two different lattice volumes V\ 
and V2, with geometries f3\ and /?2, WChPT predicts a scaling of the LECs as 



S ctt (Fi) 

(Vi) 
cf(V 2 ) 



„cff 



1 + 



1 + 



3 (A^-ft^) , {3ad 



2F 2 
4 (A 



VWV2 
V2-P2VVI) 



+ 



mc 2 S 



1 



1 

¥2 



F 2 



VWV2 



(2.17) 
(2.18) 



From the NLO partition function we can also extract information about the Aoki phase. 
The possible existence of such a phase in which flavor symmetry can be broken (with no 
analogon in the continuum theory) has been an outstanding problem for a long time. Quite 
recently in [25] the authors showed that while in the unquenched theory both scenarios (the 
Aoki and Sharpe-Singleton scenario) can be realized, the quenched theory at sufficiently 
small quark mass is always in the Aoki phase (see [40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50] 
for lattice data). 

Morally speaking in the infinite- volume case flavor symmetry breaking can occur due to the 
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fact that (Uo) ^ 0. Let us repeat here the analysis of [1] and apply it to our NLO results. 
In order to determine the value of the minimum of the potential energy we parameterize 
Uo = A + iBj ■ aj, with Gj the Pauli matrices. This makes the action in eq. (2.13) only 
depend on A. If we assume that the sign of c% is positive the potential will be a parabola, 
and the minimum is given by the parameter usually called 

" mS ™ (2.19) 



16a 2 C2 16a 2 C2 

If this parameter lies outside the range —1 to 1 then it is simple to see that the vector 
symmetry can not be spontaneously broken, and the minimum is taken by A = 1. However, 
if the minimum satisfies |e| < 1 the vacuum is determined by A* = i. As a consequence 
Bj ^ and flavor symmetry is spontaneously broken to U{1). This tells us that the region 
— 1 < i < 1 has the properties of the Aoki phase. We denote the value at which the 
transition takes place by (f 2 = jfp-- 

One can repeat the same analysis using our NLO partition function to analyze the role 
of the finite-volume corrections to this picture. The parameter e is obviously modified at 
finite V and more precisely it is given by 

mS eff 

rh ( 4 A/N a ( 9 , 9 , do\ rh ( do \ 1 



1 + —2 A(0) - — = 2a 2 d 1 + 32d 2 d 2 - 3^ - - ^ + 



16d 2 c 2 V 2F 2 y ' mVV V c 2 J a \E Ac 2 J 

The finite-volume corrections can thus change the phase of the system since they modify 
the range of the minimum of the parabola determining the potential energy of the system. 

2.3 Partition Function for Generic Number of Flavors and Fixed Index 

In this section we will consider Wilson Chiral Perturbation Theory with Nf flavors. The 
situation becomes more complicated since additional terms are allowed in the Lagrangian. 
Indeed the Wilson chiral Lagrangian for a generic number of degenerate quark flavors with 
mass m reads at LO 

Clo = ^Tr (dpUdprt) - ^Tr (u + rf) - ^Tr (u - 

WWgTt (u 2 + U^ 2 ) + a 2 W 6 Tr (u + U^Y + a 2 W 7 Tr (u - U^Y , (2.21) 
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where we have introduced a source z for the axial quark mass, and for the NLO [24] 1 it 
reads 



C 



NLO 



a Tr 
— amwQ ( Tr 



Tr 

u + u ] 



u + u ] 



+ aw^ Tr 



+a x\ ( Tr 



amwj ( Tr 



+ a X2 { Tr 



U - U ] 



amwgTr 



U 2 + U ]2 



+ a X3 Tr 
+ a 3 x 5 Tr 



E/ + E/t 



U 3 + U r > 



+ a xq Tr 



U-W 
+ a 3 x 4 Tr 
17 + 17+ 



Tr 



17 + Eft 



E/ 2 - E/ f: 



Tr 



E/-E/ f 



(2.22) 



At this point one can proceed following the same steps used in the Nf=2 analysis, namely 
expand the action up 0(e 2 ) and integrate out the pion fluctuation. Just for stylistic reasons 
we report our detailed computation in appendix B. Here we present only the final result 
and make some clarifications. 

As we can see from both the LO and NLO Lagrangians (2.21) and (2.22) a lot of terms 
appear compared to the usual Nf = 2 case. At LO there are three independent LECs, i.e. 
Wq,W7 and W$, and we recall again that in the simple case of two flavors they combine 
just in one coefficient called C2 = Wq + W%/2 while W7 doesn't enter in the computation. 
At NLO we have eleven new coefficients that are divided as follows. We have order 0(p 2 a), 
and it's written in the first line of eq. (2.22). Again one can show that for the case of 
two flavors these two terms combine in one contribution which coefficient that we called 
Co = W4 + W5/2. Then there are three terms of order 0(ma) whose coefficients are wq, wj 
and ws, and the expressions are listed in the second line of the same equation. Again for 
two flavors w-j doesn't contribute and the other terms form a combination that we called 
c 3 = w & + ws/2 previously. And finally there are the more tedious contributions of order 
C(a 3 ) with six different terms and the relative six coefficient X{ with i running from one to 
six. In the case of two flavors, only two are independent and in particular the non-trivial 
combinations d\ = x\ + £3/2 + X5/4 and d<i = x§ — 4x3 — 3xs enter in the game, whereas 
X2 and X4 don't contribute. 

Since we will work at fixed index u we define the projection to the following Fourier com- 
ponents: 



U(N f ) 



[d H U(x)\ det[E/(x) 1 



(2.23) 



Taking the Fourier sum over all components will lead back to the SU(Nf) integral. Note 
that after fixing the index the Nf = 2 case has as many terms in the Lagrangian as for 
general Nf, as there are no identities left to simplify it. At LO this partition function can 



1 Note that compared to [24] we have absorbed a factor of 2T./F 2 into our 
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be evaluated as 



'LO 



M / d H U det [U%] exp 
lu(N f ) 



mSV. 



-TV 



-a 2 VW 8 Tr 



c/ 2 + t/ t2 



a^VWe TV 



H Tr 



Uo - Ul 



a 2 VW 7 Tr 



Uo - m 



with expectation values defined by 



(F(Uo) 



1 



/LO 



j(0) 



2 LO JU(Nf) 



d H U F(U ) detpZ] e °»o 



(2.24) 



(2.25) 



using the corresponding action sjj^ from eq. (2.24). The group integral eq. (2.24) is known 
explicitly and given in appendix C. 

Since we have to take expectation values with respect to a U(Nf) integral instead of 
SU(Nf) we have derived new identities between the expectation values of the various 
terms in appendix A. 3. 

After expanding the action up to 0(e 2 ) and integrating out the fluctuations it becomes non- 
trivial to reabsorb all the terms in such a way that the NLO Lagrangian can be expressed in 
terms of the LO one. Indeed we have to make use of three relations listed in the appendix 
A. 3 to find that 



^■nlo 



Zl (m cS ,z cS ,af,af,af 



2a 3 d' 



d 2 



V daidm 



2a 3 8 2 
V da 2 - dm 



z lo (m,z,aQ,a 7 ,as) 



AX 



cff 



2 



dzdrh LO 



z lo (m,z,aQ,a 7 ,a s ) 



(2.26) 



eff -eff A eff ^eff Aeff 



a 



7 ■ 



at , as well 



where the explicit expressions for the renormalized constants m e ", z" 
as for Xf , X^, X^ s , X| are derived in the appendix B. 

When we set the source z, = the last effective coupling vanishes, X^l^o = 0. Also note 
if we were to set all the extra LECs to zero that contribute to the chiral Lagrangian at 
NLO (as it happens for SChPT in the next section), that is x\ § = = ftf6,7,8> then all 
Xf*2 3 5 = would equally vanish, and we could again write the NLO partition function 
as a LO one with the couplings renormalized through the one-loop corrections. Also the 
masses m and z would then be renormalized with the same effective LEC S cfr . This is the 
situation we find below in SChPT, and it is also true for the NLO finite- volume corrections 
in the continuum. 



2.4 Two-Point Correlation Functions for Nf = 2 

In this subsection we will calculate the two-point correlators of the scalar and pseudoscalar 
current densities, in analogy to the continuum results in [51, 52, 53]. They are defined 
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respectively by 



S (x) = ${x)i>{x) ,S h {x)=i>(x)t h i>{x) (2.27) 
P (x) = i$(x)-ys1>(x) ,P b (x) = #(x)t b75 V(x) (2.28) 

at the first non-trivial order in the e-expansion. In the previous expression tf, are propor- 
tional to the Pauli matrices t& = \ob for b = 1,2, 3. Following the standard procedure, in 
order to calculate these quantities one has to introduce the Hermitian sources s and p in 
the partition function through the replacement M — > M + sq(x) + Sb(x)tb + ipo(x) +ipb(x)tb, 
and take the following functional derivatives: 

{St{x)s < m = hsJs, M z[eM • (2 - 29) 

™ I>fyo>> - 1 s p J)SpM ZIs - p] ] '---" ' <2 ' 30) 

Before starting we recall that the isovector scalar (sb=i,2,3(x)) and the isoscalar pseudoscalar 
(po(x)) densities, as a property of the SU(2) theory, are vanishing at LO and also NLO 
and thus we will not consider these quantities. 

The additional sources lead to the following modification of the LO Lagrangian eq. (2.1) 



E, 



2 



5C LO = -s Q {x)-Tr U + U^ + i Pb (x)-Tv t b (U - W) , (2.31) 



2 



as well as to the corresponding modification of the NLO Lagrangian eq. (2.4) 

,5£nlo = s (x)ac 3 (Tr[U + t^])' - ip b (x)ac 3 Tr[t b (U - tf)]Tr[U + U^} . (2.32) 



In addition to the continuum the two-point functions are only known up to the renormal- 
ization constants Wg P , that are of the form (1 + aWg)(Sb(x)S c (0)), and likewise for the 
pseudoscalars (see [54] for details about the renormalization procedure). 
Now we have to calculate the correlators at NLO in the e-expansion. First, we expand the 
observables up to 0(e 2 ) using the NLO action, in the second step we perform the Gaussian 
integration over the fluctuations. Here we report only the results for SU (2) while we show 
explicitly the full calculations for U(Nf) at fixed index in the appendix C. 
The advantage of the fixed index averages is that we could use the compact integral rep- 
resentations for the LO partition function, eq. (C.4) or eq. (C.8) for Nj = 2, to obtain 
the NLO expressions from eq. (2.26). The logarithmic derivatives with respect to the 
corresponding couplings then generate all group averages given explicitly in appendix 2.4 
in eqs. (C.9), (C.10) from these LO and NLO partition functions. 

For practical purposes however the disadvantage of fixed index at NLO is the large number 
of LECs to enter the expressions, that is for Nf = 2 the 3 LECs Wqj^ from LO plus an 
additional 6 combinations from NLO in eq. (2.26). For this reason we have not attempted 
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to plot the NLO two-point functions calculated at fixed index. 
For the two-point function of the scalar current density we obtain 



(s cff )2 



(So(x)So(O)) 



4 



Tr[U + Ul 



NLO 



2F 2 



Tr 



Tr[<7 + Ul\ 

V \ K ' I LO 



(2.33) 



where the averages are now over constant matrices Uq £ SU(2) with the NLO or LO 
partition function eq. (2.13), respectively. Apart from the averaging partition function the 
expression in the first line completely agrees with the continuum expression, see e.g. [33], 
after setting Nf = 2 with Tr U - I) 
For the pseudoscalar sector we have 



there. 



(E eff )2 



<A(*)A(o)> 



Tr 



4 ) 

/ NLO 



T? /3 , 
+ 4M4 ( 



U + £/ n T 



+ 



2W 



Tr 



f/ 2 + C/ t2 



Tr 
4) Tr 



c/ 2 + c/ t2 



Co + Ul 



+ 4) A(x) 

/ LO 
LO ' 



(2.34) 



where we have used the su(nj) completeness relation eq. (A. 2) after summing over b = 
1,2,3. Once again the first two lines agree with the continuum expression for Nf = 2, 
apart from the different average. 

Let's start by plotting the zero-momentum correlator for different values of the mass m 
and at fixed values of C2 and lattice spacing a. More in detail we consider a hypercubic 
symmetric lattice with N L = N T = 48 (ft = 0.140461) and use F = 90 MeV and E = 250 
MeV. In table 2. we list all parameters and numerical values of the LECs used in the plots 
(for details on the explicit integrals used see appendix D). 





a(fm) 


m(Mov) 


C2(GeV 4 ) 


(il(GeV 4 ) 


d 2 (GeV 4 ) 


C 3 (GeV 4 ) 


Fig. 1. 


0.08 


7, 2, 0.5 


0.01 


10" 5 


10~ 5 


10" 5 


Fig. 2. 


0.12, 0.10 


1 


0.005 ... 0.050 


10~ 5 


10" 5 


10" 5 


Fig. 3. 


0.10 


1 


-0.1 ... 0.1 


10" 2 









Table 2. Numerical values of the lattice spacing, masses and LECs used in the 
evaluation of the two-point function from appendix D. 



In figure 1. we plot the integrated two-point correlator, eqs. (2.34) and (2.35), comparing 
our results with the continuum limit [51, 52, 53], and with the results obtained by Bar, 
Necco and Schaefer [29] (see also [30]) that are valid in the GSM* regime, where finite lattice 
spacing effects only enter at NLO, compared to LO in our LCE counting. We only display 
the pseudoscalar-correlator, at three different masses with otherwise fixed parameters: for 
the first two values chosen, m=0.5 and m=2 MeV the system is in the Aoki regime (top 
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a 3 C(t) 




Figure 1. Pscudoscalar correlators at fixed a=0.08 and different masses m=0.5 (top-left), m = 2 
(top-right) and m = 7 MeV (bottom). The different colors distinguish the continuum result, the 
results of [29] in the GSM* regime and our result. 




— 0.005 

— o.oio 

— 0.015 
— 0.020 
— 0.050 



Figure 2. Pscudoscalar correlators at fixed a=0.12 (left), a = 0.10 (right) and fixed mass m = 1 
MeV, for different values of the LEC C2 encoded by different colors. 

two plots in fig. 1.) while for the last one with m=7 MeV it is in the GSM* regime instead 
(bottom fig. 1). Our calculation agrees quite well with [29, 30] when the GSM* counting 
is valid, while it disagrees when one enters the Aoki regime. Indeed the GSM* expansion 
is then no more reliable since lattice spacing effects give LO contributions and cannot 
be considered as perturbations. Similar plots could be obtained for the scalar two-point 
function. 
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In figure 2. we plot the same pseudoscalar correlator for different values of C2 ranging from 
0.005 to 0.05 GeV 4 , at fixed m and a. In the left figure we use a lattice spacing of a = 0.12 
fm while in the right plot we have a= 0.08 fm. As one could expect for the bigger lattice 
spacing small changes of C2 bring up the parabola, making the corrections to the continuum 
more and more severe. 

Finally we want to analyze the slope of the parabola described by the correlators. In order 
to do that we can recast them into the following form given on the right hand side, 

2 l\ 
-—}N t a, (2.35) 



a 3 C(t) = J d 3 x^2(P b (x, t)P b (0)} = Ap + B P (^ 



i 



and study how the coefficient Bp depends on the LECs. As we can see clearly from figure 
2. the main correction to the continuum limit comes from the modification of the constant 
Ap that determines the value of the minimum of the parabola. Indeed when the LEC C2 
increases (we assume that at NLO the other LECs have the same effect to increase c| ), 
the minimum of the parabola becomes larger and larger compared to the continuum. 
In figure 3. we look at the value of Bp, that is the coefficient that drives the slope of the 
parabola and thus is related to the masses of the pions. Differently from what happens in 
the GSM* regime, in the Aoki regime also this parameter gets modified by lattice spacing 
effects. More in detail, plotting the value of Bp as a function of C2 it is interesting to note 
that for small enough values of C2 we are close to the continuum in the infinite-volume 
limit. Referring back to our discussion of the Aoki phase boundary in subsection 2.2, for 
— oo < C2 < c?, the system is in the so called Sharpe-Singleton scenario 2 . In that region 
flavor symmetry is not broken and all three pions remain massive. The situation changes 
when C2 is bigger than c*>, since the system enters in the Aoki phase. The value of Bp 
decreases quite rapidly to reach a lower limit of approximately 2/3 of the continuum limit. 
This is an indication that flavor symmetry is broken and as a consequence two of the three 
pions become massless. 



3 Staggered Chiral Perturbation Theory at NLO 
3.1 Introduction 

In this section we study finite-volume corrections to the LECs in the framework of SChPT. 
As we will see, the situation is easier compared to WChPT since the corrections of order 
0(am),0(a 5 ) and 0(ap 2 ) don't appear in the staggered Lagrangian. This simplifies re- 
markably the computation and allows us to write the expression for the NLO partition 
function with a generic number of flavors in term of the LO one. Let us start to re- 
view briefly some basic known facts about SChPT and its equivalence to Staggered Chiral 
Random Matrix Theory (SChRMT) in the e-regime [22]. 



Leading Order O(e ) 
Next-to-Leading Order 0(e 2 ) 
Next-to-Next-to-Leading Order 0(e 4 ) 



2 2 

m, p , a 



2 2 4 2 2 2 4 

m , mp , p , am, a p , a 



The value C*, is when ci satisfies ev = 1, see the discussion before eq. (2.20). 
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BP(C2) 




C2 



-10 



5 



5 



10 



Figure 3. The value of Bp from eq. (2.35) as function of C2. Both x- and y-axis are in units of 
10 10 MeV 4 . 



Staggered fermions are widely used to simulate quarks on the lattice. Indeed this formu- 
lation presents some clear advantages as the fact that the continuum chiral symmetry is 
not completely broken and that it is quite inexpensive to simulate numerically. However 
it doesn't solve completely the doubling problem. For every physical flavor there are four 
taste states that are degenerate in the continuum but split at finite lattice spacing because 
the taste symmetry is broken. Furthermore, in order to restore the correct number of 
degrees of freedom one has to introduce the so-called "rooting" procedure that consists of 
taking the fourth root of the fermion determinant. 

The effective chiral Lagrangian that describes the staggered formulation including finite 
lattice size corrections has been introduced in [6] for the one flavor case and generalized to 
multiple flavors in [7]. The authors added all 0(a 2 ) terms to the continuum Lagrangian 
that are compatible with the staggered symmetries. As usual the breaking of the chiral 
symmetry from G = SU(4Nf)n x SU(4Nf)L to SU{ANj)y is associated with the existence 
of light (pseudo) Goldstone boson fields that are collected into a ANf x ANf unitary matrix 



U. For example in the case Nf = 3 the matrix U £ SU(12) can be parameterize as 

/ U 7T+ K + \ 



where = u, tt + , K + ... are the 4x4 matrices that take into account the taste degrees 
of freedom and that can be written in the Dirac basis as U = UbT^, with denoting 
Tb = {£5, ^£5, i£ M £„, £/}. In the LCE regime the LO Lagrangian, which is of order 
(D(p ,m,a ), reads as [7] 



Table 3. Contributions to the Staggered Chiral Lagrangian in the Aoki regime. 
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Clo = ^Tr (d^Ud^ - |Tr (m*U + U^M) - a 2 CiTr (^75^75 



4 



—a 



J [Tr (Efy) Tr (t/^)] - a 2 ^ £ [Tr (tf 7/lB ) Tr (fft 7/lB 



-a 2 C 6 ^Tr^^V) ; (3.!) 

where as usual F and S are the pion decay constant and the chiral condensate, respectively, 
while the 4iVy x 4Nf matrices 7^ are the generalizations of the ordinary 4x4 Dirac ma- 
trices £^ (see [7] for details) 3 . In addition to the continuum Gasser-Leutwyler Lagrangian 
there are some taste-breaking contributions and as a consequence some new LECs usually 
denoted as C±, C2A1 C2 V> C3, C4, C^a, C5 V, Cq. For the one-flavor case the situation sim- 
plifies since all the two-trace terms in eq. (3.1) can be Fierz transformed into one-trace 
terms (see [6] for details). As we can see from table 3., if we want to go beyond LO, at 
NLO we will only get effects from one-loop finite- volume corrections of order C(e 2 ). Fur- 
ther LECs arising from the discretization effects of order 0(a 4 ), 0(a 2 p 2 ) and 0(a 2 m) will 
only appear at NNLO, being of order 0(e 4 ). The corresponding terms are carefully listed 
in [55]. 

Quite recently it has been shown in [22] that in the e-regime SChPT is equivalent to 
SChRMT, including all one- and two-trace terms. In the latter theory the taste braking 
terms are introduced by adding a taste diagonal matrix to the usual Dirac operator 

V=( .L iW )®h + a 2 T , (3.2) 



1 







with W a random matrix of size (N + v) x N with complex entries, and where the explicit 
form of T and its relations with all the terms in the taste breaking potential are given 
in table 1. of [22]. This correspondence is nothing else that the analogous relation for 
the staggered case between the e-regime of WChPT and WChRMT analyzed previously, 
valid at LO in the e-expansion. The aim of this section is to study the finite-volume 
correction to the LECs Cj following the same procedure utilized in the previous section 2 
for WChPT. As a byproduct of the calculation we will find how the taste splittings, i.e. the 
difference between the mass square of the non-Goldstone boson and the Goldstone bosons, 
are modified. 

3.2 Staggered Chiral Partition Function at 0(e 2 ) for General Nf 

In order to calculate the partition function at 0{e 2 ) in the staggered case we follow the 
same steps used in the analysis of the Wilson chiral Lagrangian. Let us start to rewriting 



There is also a mass term for the taste singlet pion that we have dropped. 
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the partition function as 



d H U(x) e~ s = / D H U e~ Su o Z^(U ) , 

SU(4N f ) JSU(ANf) 



(3.3) 



where we have divided as usual the integration over the zero-modes Uq from the integration 
over the fluctuations £ and where 

Expanding the function Z^{Uq) up to order 0(e 2 ) one can perform the Gaussian integrals 
over the fluctuations (for details about the expansion and the integration see appendix 
E). The next step is to reabsorb the finite- volume corrections into the LO Lagrangian by 
re-exponentiating all the terms found in the appendix E. At the end we can conclude that 
the NLO order partition function can be rewritten as the LO partition function with some 
renormalized LECs 

(3.5) 



^NLO = -Jj ^LO y-> , t>i 

These renormalized LECs following from the calculation in appendix E are shown in table 4. 
Similarly to the Wilson case the renormalized effective constants depend on the geometry 
of the system through the propagator A(0). As a further consequence we have extended 
the equivalence between SChPT and SChRMT up to NLO, as the form of the LO partition 
function is preserved for any number of flavors. 



AF 2 N 



A(0)J 



Cfy — Civ 



C 2V (16N 2 -2)+4C 3 N f 



2N f F' 2 



A(0) 



ff _ C 3 (16^-2)+2C 2 yJV f 
^3 - °3 2N f F 2 A W 



C§ = C 5V (l - ^A(O) 



^A(O) 



cf = C fi [ 1 



cf 



Ci 1 



fA(0) 



ca=^- c "" 8 S:r C4K/ A(o) 



/~<cS /"i 

4 — O4 



2N f F 2 



C i (l6N 2 -2)+2C 2A N f 
2N7W A (°) 



C? A = CfiA ( 1 



8AT 



fA(0) 



Table 4. The renormalized LECs in SChPT. 



From the previous computation one can immediately understand how and if the finite- 
volume corrections affect the taste symmetry. Usually to study the taste symmetry viola- 
tion one looks at the taste splitting A in the pion sector, i.e. the difference between the 
mass square of a non-Goldstone pion and of the Goldstone one. At LO this quantity can 
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be derived from a tree level expansion of the chiral Lagrangian, and indeed the masses of 
the non-neutral meson 4 composed of quark b and c can be written as 

m 2 = fi(mi, + m c ) + a 2 A^ B . (3-6) 

The terms A are related to the LECs through the relations [7] 

A P = , (3.7) 

1 6 

Aa = (Ci + 3C 3 + C 4 + 3C 6 ) , (3.8) 

1 C 

A T = ^ (2C 3 + 2C 4 + 4C 6 ) , (3.9) 
16 

Av = -^ (Ci + C 3 + 3C 4 + 3C 6 ) , (3.10) 

Aj = H (4C 3 + 4C 4 ) . (3.11) 

These splittings concern the Pseudoscalar (P), Axial- Vector (^4), Tensor (T), Vector (y) 
and Singlet (/) taste pions respectively. Since the LECs are modified at finite-volume the 
taste splitting get modified as follows 

A^ LO = A P = 0, (3.12) 

A n,o = Aa _ « J 8Ar/[Ci + 3C6] + lC i + 3 C :mNr ^C 2v+ C 2Al N^ m 

(3.13) 



Ag"> = L 2NfC , + [C 3 + C 4l(1 ^-^ +2 ^ + C, <W A(0) (314) 



A ™ _ Av _ w 8iV/[Ci + 3C6] + ^i^-yfa A(0) 



(3.15) 



A NLO . Aj _ |2 / [C 3 + GJQflV? - 2) + 2 [Qy + C2 ,]ivA a(q) (s i6) 



4 For flavor neutral mesons the situation is more complicated and other terms have to be introduced in 
the chiral lagrangian. 
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4 Summary and Discussion 



In this paper we have computed 0(e 2 ) finite-volume corrections in the so-called e-regime 
as they arise in Wilson and Staggered Chiral Perturbation theory. Thus we have taken 
into account both corrections to infinite-volume and to the continuum limit. In SChPT 
0(a 2 ) effects only enter at LO, parameterized through a large number of in total 6 new 
low-energy constants. In contrast in WChPT such 0(a 2 ) effects enter both at LO and at 
NLO in the e-expansion, leading to a total of 3 plus 9 LECs, respectively. 
In consequence SChPT, although more complicated at LO, will remain simpler at NLO. 
In particular on the level of partition function the effect of NLO can be entirely expressed 
by renormalizing the LO LECs with corrections that we explicitly computed. As a second 
consequence the known equivalence between SChRMT and SChPT continues to hold at 
NLO. The drawback of the staggered formulation however remains, that the corresponding 
SChRMT has not been solved analytically to date. In addition our results for the NLO 
LECs provide us with the finite- volume corrections to the taste splittings as an application. 
Turning to the Wilson case much is known about the LO spectrum of the Wilson Dirac 
operator at fixed index, due to the equivalent WChRMT picture. This equivalence breaks 
down at NLO for an arbitrary number of flavors at fixed index including Nf = 2, as ex- 
tra derivatives appear when trying to express the NLO partition function through the LO 
partition function with renormalized, effective couplings. Only for the original zero- mode 
group integral without fixing the index, and for the special case of SU(2), the NLO partition 
functions keeps its functional form. As a consequence we can use the two corresponding 
effective couplings to quantify finite-volume and finite-lattice spacing effects on reaching 
the Aoki-phase. For the same reason the finite- volume and 0(a 2 ) corrections remain sim- 
plest for the two-point functions for SU{2), which we have computed in the scalar and 
pseudoscalar case explicitly, and plotted in the pseudo-scalar sector for illustration. 
Let us comment on finite- volume corrections to the positivity constraints on individual and 
certain combinations of LECs. At LO these were based on the positivity of the partition 
function at fixed index [11], on Hermiticity arguments for the generating functional for 
the spectral density of Wilson Dirac eigenvalues [25], and on the mass split using partially 
quenched WChPT [23]. It appears that neither line of argument can be easily translated to 
NLO, by simply replacing the LECs by effective ones. This has to do with the fact that at 
fixed index and/or for Nf > 2 the functional form of the NLO partition function changes 
compared to LO. 

In principle, the effect of NLO on the Wilson Dirac spectrum could be computed in the stan- 
dard way, by introducing graded or replicated partition functions as generating functionals. 
However, due to the loss of determinantal structure of the partition function at NLO, that 
is observed at fixed index and certain vanishing couplings at LO (see appendix 2.4 where 
we also computed an extended version), and due to the loss of the WChPT- WChRMT re- 
lation, this seems to be a formidable task. Such a result would be very interesting in order 
to explain asymmetric effects on the spectrum attributed to NLO corrections in [17, 19]. 
On the other hand the extension of our computations to vector and axial-vector two-point 
correlation functions should be straightforward, both for SU(Nf) and for U(Nj) at fixed 
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index. The practical difficulty for the resulting expressions (except for SU(2)) however is 
then to determine the large number of effective LECs from actual data, on which these 
quantities will depend. 
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A Zero Mode Group Integral Identities 



A.l General SU(N f ) Case 

In order to derive some SU (Nf) group identities among the expectation values of the various 
trace terms we follow the strategies adopted in [33, 51]. We introduce the differentiation 
with respect to the group elements U^i of SU(Nf) defined as 



V 6 = i(t b U) 







kl 



dU kl ' 



(A.l) 



where t b are the generators of the algebra su(Nf) and satisfy the completeness relation 

{tb)ij{tb)ki = 7^ y>u8jk ~ j^SijSkij . (A.2) 

This leads to the following derivatives 

V b U = it b U , V 6 i7 t = -iUh b . (A.3) 

Considering that the Haar measure is left invariant, the integrals over total derivatives with 
respect to V& have to vanish and thus for example 







SU(N f ) 



d H U V c { Tr[t c G{U)] exp ( ^-Ti 



U + U ] -a 2 W 6 VTr( U + U ] 



a 2 W 7 VTr 



U-W 



a 2 W 8 VTr 



U 2 + W 2 



(A.4) 



holds for any choice of the function G(U). Throughout this appendix U = Uq is a constant 
matrix and for simplicity we drop the subscript compared to the main text. The following 
brackets denote the expectation value with respect to the integrand 

(F(U)) = (A.5) 

f dnU F (jjy^^[U+U^}-a 2 W a VTr([U+U^]) 2 -a 2 W^ 

JSU(Nf) 
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Choosing G(U) = U — we obtain the following identity 
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-2a 2 W 6 V (Tr U + U ] 



U 2 - U r ' 



-2a 2 W 7 V ( (Tr 
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Tr 
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J N f \ L J/ L 
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U 2 - U r ' 



Tr 



17- [7 f 



(A.6) 



A.2 The 5C/(AT / = 2) Case 

Here we rewrite the previous identity for the particular and more simple case of Nf = 2 
quarks. Indeed in this case some relations between the trace terms can be used to simplified 
considerably what we have found in (A.6). More in detail for any matrix U that belongs 
to the group SU(2) the following relations are valid 

Tr[U - [7 f ] = Tr[U 2 + U^ 2 ] = - (Tr[U + [7 f ]) 2 - 4 , 



(A.7) 



Tr[U 3 + C/t 3 ] = I (Tr[U + - 3Tr[U + W] 



Note that these identities hold without taking an expectation value. Thus one can see that 
the general Nf expression of the identity reduces for this case to 







- + 16a 2 c 2 V ) (Tr 



U + C/ f 
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(Tr 



U + C/ f 



a 2 c 2 V( Tr 



U + U ] 



-4mSF (A.8) 
that is the analogon of eq. (A.6), and we recall that we denoted by c 2 = Wq + W%/2. 
A. 3 General U(N f ) Case 

In this part of the appendix we derive U(Nf) group identities that will be useful when 
working at fixed topology. Differentiation with respect to the group elements U^i of U £ 
U(Nf) is defined as before 



V fe = i(t b U) 







ki 



dU k i ' 



(A.9) 



where tj, are now the generators of the algebra u(Nf) and satisfy the completeness relation 



{tb)ij(tb)kl = ^5il&jk- 



(A.IO) 



Because of changing from SU(Nf) to U(Nf) group integrals we will have an extra factor 
det[C/] y included in the integrand, the derivative of which reads 



V{,det[£/] = iTr[t 6 ]det[l7] . 



(A.ll) 
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Again the integrals over total derivatives with respect to V& has to vanish: 

zT,V 
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(A.12) 



for any choice of G(U). Here we have added an extra source term z for later convenience. 
The brackets denoting expectation values are now labeled by the index v in order to 
distinguish from the previous subsection. Also we use the following abbreviations: 



mSV = m, zT,V = z, aj = a 2 WjV for j = 6, 7, 8. 



We now derive a series of identities. Consider 
• t c = t 5 cfi and G[U] = 1: 
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(A.14) 



summing over i c and G[U] = U — W: 
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(A.15) 



We found two further identities with up to cubic powers of U and by choosing G[U) = 
U + . However, these identities contain new terms not present in the equations we wish 
to simplify and hence they are not useful. 



B Wilson Chiral Perturbation Theory for General Nf at Fixed Index 

In this section we focus on fixed index, the reason being that we then can compute the 
group integrals more explicitly. Also we have more group integral identities available from 
the previous subsection. In this way we can express the NLO partition function through 
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the LO one at renormalized couplings and derivatives thereof. 

The partition function up to 0(e 2 ) can be written as a sum of two contribution and 
that read 
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They are the contributions of order O(e ), and 0(e 2 ) respectively, where we have split the 
former one into the zero-mode and propagating mode part. The order 0(e) vanishes. 
Now we can proceed following the same steps used in the analysis of the two-flavor theory. 
We begin rewriting the partition function for the Nf flavors with fixed index v as 



U(N f ) 
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The function Z^(Uq) can be calculated expanding eq. (B.4) up to 0(e 2 ) order and then 
performing the integral over the Gaussian fluctuation £(x). One obtains 
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Using the relations (A.13)-(A.15) from the previous appendix A. 3 we can rearrange several 
terms present in eq. (B.5) as a sum of other contributions. After several manipulation the 
final answer reads: 



— A(0) H x$Nf + x 4 -x 5 — - ) + a 3 V x 5 + x 6 H 

AT f 2F 2 V ; 2W 8 V W 8 J V W 8 



ig - 1 zVT. 
N f 2F 2 



av ( ( W 7 



16aW 2 



Tr [U - Ul 



( W8 ^lzl +W6+Wr ) A(0)+ammV a ^ m] 



F 2 V"" 
2a2 V („ r N f ~ 1 N fWs - 2W 7 \ 

— r— + ^ j 

2a 2 V / , A"? - 1 NtWg - 2W 6 



4W 8 



Tr[c/ 2 + C/ t2 ] 



A(0) +amu) 6 V 



F 2 



W 7 - 



N 



f 



2N 



f 



A(0) + amV { w 7 - I x 4 - x 5 —) -^—) 
5 W a J 4W 8 J 



W 6 + W 7 



a 6 V 23 



W 8 
X5W6 
Ws 

(W 6 + 2W 7 ) \ azY.V 



W 7 
'Wg 



) (Tr \u -Ul\) Tr ([/ + u\ ) 



Tr (t/ 2 + Ul 2 ) Tr ([/ + U 1 ) - a 3 xi V Tr (u + ) ' 
Tr (t/ - t/J ) Tr (f/ + J/t ' 



16aW 2 



Tr [u + Ul] 



[Tr [U - Ul 



x 4 — XK- 

W 8 J 4W 8 

At this point it turns out to be useful to define the new renormalized masses and LECs as 
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Note that at NLO the quark and axial quark masses do not renormalize with the same 
LEC S any more. The renormalized constants m , m , tigy 8 contain both 0(1) and 0(e 2 ) 
parts, the constants Xf^z ar e all 0(1) only, and Xf is of 0(e 2 ) only. This is because of 
the following NLO expression for the partition function: 
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Scalar and Pseudoscalar Currents in WChPT at Fixed Index 



(B.16) 



In this appendix we will complement the main body of this paper by computing the par- 
tition function and scalar and pseudoscalar two-point functions for an arbitrary number 
of flavors Nf at fixed index v. We begin with the partition function. Given the previous 
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appendix we only need to compute it to LO as the NLO one can be expressed through it. 
It is defined as 
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with the rescaled quantities rh = mVT,, z = zVE, a 2 VWj = d| for j = 6,7,8. We have 
dropped the index of Uq here and in the following. This integral has been calculated in 
the literature in a series of works [10, 11]. Let us briefly review and slightly extend their 
results. 

Consider the following group integral that contains the above case for clq = a-j = 0: 



jN f ,u = 



d H U de^C/fexp 



£(a/I¥[^]+/3/IV[[/ti] 



(C.2) 



After diagonalizing the matrix, U — > v diag(e* 6 ' 1 , . . . ,e ldN t)v^ , with v € U(N f) /U(l) N f , it 
can be written as a determinant over a single integral: 



N f 
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In the first step have rewritten the absolute value square of the Vandermonde determinant, 
the Jacobian resulting from the diagonalization, and pulled the exponential prefactors into 
the respective determinants. In the second step we have applied one of the de Bruijn 
integration formulas. The constant Cjv* is the volume of the coset integral over v. 
As a consequence for d& = dj = we can write the corresponding A^-flavor partition 
function as a determinant of a single flavor partition function, 
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where 



Z^q l ' u (rh, z, 0, 0, dg) = C\ I exp[i9u + mcos (9) + izsin (9) — 2d\ cos (26*)] . (C.5) 
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In the particular case of two degenerate flavors eq. (C.4) reduces to 



Z N t =2 >" ~ (Z N f =1 ' v f - z N f =1 ^ +1 Z N f 
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(C.6) 



where we have suppressed the arguments. The remaining two LECs clqj can be switched 
on by performing two Gaussian integrals on the above formulas, following [11]. In addition 
it has been shown in [25] based on Hermiticity that both LECs have to be non-positive, 
06i07 < 0. Consequently we obtain the following expression: 
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.oo 167r|a 6 a 7 | 

Now let's derive the partition function at fixed topology in the case of two flavors. In fact 
the Gaussian integration over the partition function (2.22) can be performed explicitly, and 
we obtain that 

Z^o(rh,z,a 6 ,a 7 ,a 8 ) = y J* d0 x dB 2 e i( - 9l+62 > (l - e^ 1 "*")) e A ( cose i+ cose 2)+^( sine i+ sin ^) 

x e 4S s ( cos Sl + cos ^ 2 ) 2 e -4 "? ( sin ei + sin ^ 2 ) 2 e — 2 "l ( cos + cos 2 ^ 2 ) (C 8) 

Although in this expression the integrals do not factorize it is very useful for a numerical 
integration. The normalization constant C2 is not important as it drops out in expectation 
values. 

We can now compute the NLO two-point functions for general Nf at fixed index, following 
the same lines as in the previous appendix B. In fact written in terms of group averages 
(. . .)", where the superscript denotes the index, the expressions with or without fixing the 
index don't differ. This is because the propagating modes that we contract always live in 
SU(Nf). The only difference is that for fixed index at Nf = 2 we no longer have the SU (2) 
identities at hand, e.g. Tr[U — W) / no longer applies. 

For simplicity we will only present that unflavored scalar and flavored pseudoscalar two- 
point functions as in the main text. The results we obtain are 
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Note that in deriving the expression for the zero-momentum correlation functions, which 
are functions only of the Euclidean time t, we would have to make use of the relations 
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D Explicit Computation of Partition Function and Currents for SU(2) 

In this appendix we will derive explicit integral representations of the scalar and pseu- 
doscalar current densities whose formal expressions are given in subsection 2.4. Since we 
are dealing with the two-flavor case we can describe the group manifold using the familiar 
parameterization of SU (2) 

Uq = (cos a + in ■ a sin a) (D.l) 

where h is a three-dimensional unit vector, the u's are the Pauli matrices and < a < 
2tt. With this parameterization for an arbitrary element of SU(2), the normalized group 
measure is 

J dU = J dn n J dasin(a) 2 . (D.2) 

The partition function can thus be written in a more manageable way using this parame- 
terization as 

^nlo = — 9 / da sin (a) 2 exp 2mT, eS V cos (a) — 16a 2 c| ff V cos (a) 2 , (D-3) 
2ir z J L J 

and correspondingly for LO by dropping the superscript eff and having a different normal- 
ization constant C. As a check we obtain for C2 = the known result for the equal mass 
SU{2) partition function Z^o = CI\ (2m) / (2rhir) in terms of a modified Bessel function. 
The expressions for the two-point scalar and pseudoscalar current correlators derived in 
subsection 2.4 can be rewritten as 
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E Staggered Chiral Perturbation Theory for General Nf 



In the following we report explicitly all the terms arising in the e-expansion up to order 
0(e 2 ) of the partition function Z^(U ) defined in the SChPT subsection 3.2. At LO O(e ) 
one obtains 
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while for the first order O(e) we have = , and for the second order one finds 
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Now one can perform the Gaussian integrals over the fluctuations, and one finds that 
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g^p2 A(0)y ^ [Tr (C/oTms) Tr (C/ 75 M ) + ^- c -] 

A(o)y E M (c/ 0>5 ) Tr (c/t 75 ' 



Because in this appendix we didn't have to use any group integral identities the same 
relations hold for fixed index, by adding detft/g] inside the zero- mode group integral. 
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